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ABSTRACT 

A  general  study  of  conductor  backed  coplanar  waveguide  is  presented.  The 
impedance  (Z0)  and  effective  dielectric  constant  (z„ff)  of  conductor-backed  coplanar 
waveguide  (CBCPW)  have  been  calculated  by  using  a  variational  method  and  the 
boundary  point  matching  method.  In  this  present  work  only  the  TEM  dominant  low 
frequency  propagation  mode  of  coplanar  line  has  been  considered.  Experimental  facili- 
ties are  vector  network  analyzer  (HP8409)  and  bench-instrument  measurements. 
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I.     INTRODUCTION 

A.  BACKGROUND 

In  current  microwave  integrated  circuit  design  practice,  microstrip  transmission  line 
is  almost  universally  employed  for  circuit  construction.  Although  convenient  and  easily 
fabricated,  microstrip  line  has  the  disadvantage  of  relatively  high  ambient  field  and 
coupling  to  outside  circuits.  Further,  the  incorporation  of  semiconductor  devices  and 
lumped  elements  into  microstrip  circuitry  is  difficult  in  that  a  via  hole  through  the 
substrate  must  be  introduced,  in  order  to  effect  line  to  ground  connection  of  the  device. 

Coplanar  transmission  line  has  the  advantage  that  both  line  and  ground  conductors 
are  available  on  the  upper  surface  of  the  circuit,  and  devices  are  easily  incorporated.  In 
addition,  the  spurious  coupling  between  coplanar  lines  and  adjacent  circuitry  is  known 
to  be  significantly  lower  than  with  microstrip.  It  is  known,  however,  that  spurious 
propagation  modes  and  anomalous  losses  may  occur  on  coplanar  lines  as  a  result  of 
improper  mode-excitation  of  the  line,  often  due  to  incorrect  termination  of  the  coplanar 
ground  planes  [Ref.  1,2,3], 

B.  COPLANAR  LINE  TERMINOLOGY 

The  term  coplanar  line  is  generally  applied  to  planar  transmission  systems  in  which 
a  ground  conductor  accompanies  the  line  that  is  active  relative  to  ground  return,  on  the 
surface  of  the  dielectric  substrate.  In  its  original  definition  by  Wen  [Ref.  A],  the  term 
coplanar  line  referred  to  a  system  composed  of  a  planar  transmission  line  accompanied 
on  two  sides  by  ground  planes  (Figure  la).  The  dielectric  substrate  was  assumed  to 
have  effectively  infinite  thickness.  The  term  coplanar  strips  is  sometimes  applied  to  a 
two-conductor  system  of  line  with  single  ground  return  (Figure  lb). 

The  widespread  use  of  microstrip  line  for  microwave  integrated  circuit  construction, 
with  a  ground  plane  backing  on  the  dielectric  substrate  (Figure  2  on  page  2a),  has  led 
to  interest  in  conductor  backed  coplanar  (CBCP)  line,  which  can  be  integrated 
compatibly  on  the  same  substrate  with  microstrip  (Figure  2  on  page  2b).  The  CBCP 
line  is  also  sometimes  called  conductor  backed  coplanar  waveguide  (CBCPW). 

The  ground  plane  of  CBCPW  adds  strength  to  a  semiconductor  chip  used  in  micro- 
wave integrated  circuit  construction  and  provides  a  low  loss  ground  return  at  microwave 
frequencies.  The  conductor  backed  coplanar  line  is  the  principal  object  of  the  present 
research.     To  date,  only  very  limited  literature  reports  have  been  published  on  the 


transmission  line  parameters  (characteristic  impedance  Z0  and  effective  dielectric  con- 
stant E„ff)  [Ref.  5,6]. 


Copper  plane 


Substrate 


Figure  1. 


Coplanar  systems 


Copper  plane 


Back-copper  plane 


Figure  2.      Conductor  backed  planar  waveguides 


In  the  present  work  only  the  TEM  dominant  low  frequency  propagation  mode  of 
coplanar  line  will  be  considered.  The  Transverse  Electric  (TE)  and  Transverse  Magnetic 
(TM)  hybrid  modes  which  can  propagate  on  planar  systems  will  not  be  considered.  The 
restriction  to  the  assumed  TEM  mode  has  been  employed  in  numerous  investigations 


of  microstrip  and  other  planar  systems,  and  has  been  found  to  provide  useful  trans- 
mission parameters  at  frequencies  up  to  10  GHz. 

An  additional  characteristic  of  propagation  on  multiconductor  transmission  systems 
is  that  odd  and  even  modal  excitations  may  propagate,  resulting  in  undesired  terminal 
characteristics  for  the  line.  The  current  flow  characteristics  of  the  dominant  TEM  mode, 
and  a  principal  even  mode  are  shown  in  Figure  3.  Suppression  of  the  unwanted  modes 
is  normally  possible  through  adequate  ground  terminations  for  the  coplanar  ground 
planes,  and  care  with  the  input  excitation  of  the  line. 


TEM     MODE 


TOTALLY     EVEN     MODE 


Figure  3. 


Current  flow  characteristics  of  propagation  modes 


II.     ANALYSIS  OF  TEM  PROPAGATION  MODE 

A.     REVIEW  OF  ROWE  &  LAO  BOUNDARY  POINT  MATCHING  METHOD 

In  this  section,  the  boundary  point  matching  method  for  a  CBCPW,  which  was 
presented  by  Rowe  and  Lao  [Ref.  7],  is  reviewed,  and  the  impedance  (Z0)  and  the  effec- 
tive dielectric  constant  {zreff)  for  the  conductor  backed  coplanar  line  on  a  dielectric 
substrate  are  calculated.  The  coordinates  employed  in  the  analysis  are  given  in 
Figure  4  on  page  5.  The  line  is  assumed  to  be  enclosed  in  a  shielded  enclosure  of  width 
21,  open  at  the  top.  The  enclosure  width  will  be  taken  to  be  large  enough  to  cause  neg- 
ligible perturbation  of  the  line  fields. 

To  find  the  solution  of  Laplace's  equation,  define  the  two  regions  shown  in 
Figure  4  on  page  5: 

1.  The  lower  region  is  x,  y  coordinate  system  (in  the  dielectric). 

2.  The  upper  region  is  x',  y'  coordinate  system  (in  air). 

The  solution  of  Laplace's  equation  in  the  upper  region  is: 


0  Vy)  =  V«;  cos(k,x')  exp(  -ty)  (1) 


where. 


A-(2/-1)g  (2) 

2/  y  } 


i  =  Y  +  D  W 


Assuming  width  D  equal  to  C/2  produces  negligible  effect  on  the  calculated  parameters. 
The  solution  of  Laplace's  equation  in  the  lower  region  is: 

oo 

4>~(xy)  =  2_J>i  cos(A>v)  sinh(A^)  (4) 

The  boundarv  conditions  at  the  surface  of  the  dielectric  are: 
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Figure  4.      Geometry  of  a  conductor  backed  CPVV 
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4>+(0,0)  =  Yj°i 


/=1 


0-(O,//j  =  £^sinh(A///) 
0+(O,O)  =  010,//) 


(5) 

(6) 

(7) 


therefore, 


at  =  bt  sinh(A///) 


(8) 


A  B 

In  the  gap  (  —  <  x  <  —  )  between  the  conductors,  and  outside  of  the  ground  con- 

C 

ductor  (  —  <x  <  I ),  the  surface  normal  component  of  the  displacement  field  is  contin- 


uous. 


Dt  =  D 


n        "n 


(9) 


where, 


Dn  =  e 


C(f) 

oy 


(10) 


Apply  the  boundary  condition  at  the  surface  of  the  substrate  to  equations  (9)  and  (10). 


8(f)+ 


y'=0         —         ErE0 


6<b~ 


CV 


y=H  ' 


bi  =  ~ 


sinh(A(//) 


(11) 


therefore, 


2^a(-  cos(A'^){l  +  Ercoih(kiH)}  =  0, 


4<.<- 


B 


C 


<x<  I 


(12) 


where  cr  =  relative  dielectric  constant  of  substrate. 

Generally  the  voltage  on  the  inner  conductor  is  set  to  unity  so  that  the  conductor 
capacitance  is  numerically  equal  to  the  charge  on  this  conductor.  The  unity  potential 
on  the  center  conductor  and  the  zero  potential  on  the  outer  conductors  lead  to 


oo 

2_at  cos(  lip) 


.0, 


^  A 

B    .      ,  C 


(13) 


The  above  Fourier  summations  are  terminated  at  i=N  to  make  problem  suitable  for 
numerical  analysis.  To  find  the  N  unknown  a,'s,  x  takes  on  N  discrete  values  on  the 
dielectric-air  interface;  at  coordinates: 


*;  =  7rn~>       y=l,2,....,iV  (14) 

Equation  (12)  and  (13)  for  the  a,'s  can  then  be  written  as  an  N  by  N  matrix  equation. 

N 

^jnjfl^dj,        j=\,2,...,N  (15) 


/=i 


where, 


A  B  C 

Ckj  cos{kpcj){  1  +  crcoth(/c^)},        —  <  Xj  <  —  ,      —<Xj<l 

mji=\       ft    >  **  *  ,     ,  c  (16) 

(.co s(ktxj) ,  Xj  <  —  ,       —  <  Xj  <  — 


and 


fl.        -V;<4 

lO,  Xj  >  f 

The  matrix  equation  (15)  is  solved  by  Gaussian  elimination  to  find  the  a/s. 

In  order  to  compute  the  impedance,  the  capacitance  per  unit  length  of  guide  from 
the  center  conductor  to  ground  is  calculated  by  integrating  the  normal  component  of 
D-field  on  a  Gaussian  surface  enclosing  the  center  conductor  only.  Figure  5  on  page  8 
shows  the  surface  S  chosen  for  the  integration.   The  capacitance  is: 


C  =  £0    ErEdS  =  —  £0    ErV(f)dS 


(18) 


therefore, 


C=-2c0^a/sin(/r/-^-j^-){l  +  crcoih(k,H}}  (19) 


1Q-9 

36n 


where  £0  =  4r~  =  8.84  x  10"12 


Figure  5.      Gaussian  surface  in  the  capacitance  calculation 

The  static  impedance  and  effective  dielectric  constant  can  be  calculated  from  stand- 
ard TEM-modal  relations: 


^o  = 


1 


VyJ  CCQ 


(20) 


and 


C 


(21) 


where  v  =  3  *  108  m/s,  C  is  the  capacitance  with  the  dielectric  layer  present,  and  C0  is 
the  capacitance  when  the  dielectric  is  replaced  by  air. 

Figure  6  on  page  9  and  Figure  7  on  page  10  show  the  calculated  results  of  impedance 
and  effective  dielectric  constant  of  CBCPW  using  boundary  point  matching  method. 


ROWE  AND  LAO  METHOD 
IMPEDANCE 
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Figure  6.       Calculated   characteristic   impedance   of  CBCPW   by   boundary   point 
matching  method:     (c,  =  3.52.  H  =  1.5  mm) 


ROWE  AND  LAO  METHOD 
EFFECTIVE  DIELECTRIC  CONSTANT 


10 


A/H 


Figure  7. 


Calculated  effective  dielectric  constant  of  CBCPW  by  boundary  point 

(e,  =  3.52,  H  =  1.5  mm) 


matching  method: 
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B.     REVIEW  OF  YAMASHITA'S  VARIATIONAL  METHOD 

A  method  for  computing  impedance  (Z0)  and  effective  dielectric  constant  (cre//)  of 
conductor  backed  coplanar  waveguide  is  based  on  the  application  of  Fourier  transform 
and  variational  techniques  [Ref.  S].  Figure  8  shows  the  geometry  of  CBCPW  for  vari- 
ational method. 


AirU,  =  1    ) 
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A 

Figure  8.      Geometry  of  CBCPW 


The  static  potential  distribution  (f>(x,y)  in  the  planar  line  structure  satisfies  Poisson's 
equation: 


v^(.^)  =  -|pM 


(22) 


where  p(xy)  is  charge  distribution  on  the  surface  of  the  conducting  strips. 
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A  variational  expression  for  the  line  capacitance  is: 


p{xy)4>{xy)dl 
where  Q  is  the  total  charge  per  unit  length  of  strip: 


Q2 

C  =  - (23) 


Q=\p{x)dx  (24) 


We  assume  that  the  strip  is  infinitely  thin.    For  the  infinitely  thin  strip  case,  the  charge 
distribution  assumes  the  form 

p(xy)=J[x)6(y-H)  (25) 

where  t{x)  is  even  function  and  S(y  —  H)  is: 

S(y  -  H)  -  |0j       y¥=H       .    Dirac,s  5  function.  (2°) 

Now  define  the  Fourier  transform  4>(fiy)  of  4>{xy)  ■ 

<KM  =        4>{xj)e!pxdx  (27) 

J-oo 

With  use  of  the  known  Green's  function  for  this  problem  [Ref.    9  ],  the  Fourier  trans- 
form of  the  potential  distribution  at  y  =  H  is: 

w^=eoimi+£,coth(i/M//)}  (28) 

2nd1 
C  =  — " (29) 

f(P)4>{fl,H)df3 

<o 

Insert  the  equation  (28)  into  (29),  then  the  final  variational  expression  is: 

12 


c  = 


™nQ 


■o^. 


f 


(1  +  ercoth(pii)}0 


(30) 


dfi 


where. 


f{H)  =  \Ax)eipxdx 


(31) 


where  f(x)  is  the  density  distribution  of  the  charge  on  the  strips,  and  symmetry  with  re- 
spect to  x  =  0  is  assumed. 

As  shown  in  Figure  9  and  Figure  10,  two  forms  of  elementary  charge  distribution 
on  the  center  conductor,  f(x)=  constant  and  l\x)=  |x|,  were  tested.  The  negative  charge 
on  the  two  ground-plane  conductors  is  assumed  to  be  constant  in  both  cases.  The 
computation  was  found  to  be  insensitive  to  the  form  assumed  for  the  charge  distribution 
on  the  top  ground  planes,  hence  the  uniform  (negative)  charge  distribution  was  assumed. 


Figure  9.      Assumed  step-charge  density  function  of  CBCPVV:     (  f(x)  =  contant  ) 
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In  case  of  Figure  9  on  page  13,  the  assumed  elementary  charge  distribution  function  is: 


I,      — j-s*sT 
Ax)-|_Ml      fsws-f 

0,  otherwise 


(32) 


The  Fourier  transform  of  f(x)  is: 


/w-f 


sin(/i  -4- )  +  iff  MP  -j )  -  sin(j8  -£- ; 


} 


(33) 


and, 


Q  =  A  +  M{D  -  C) 


(34) 


Figure  10.      Assumed  linear-charge  density  function  of  CBCPW:      (f(x)  =  |x| 


In  case  of  Figure  10,  the  assumed  elementary  charge  distribution  function  is: 
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/w  = 


w, 

A     .      .  A 
2              2 

-M, 

-<k|<T 

o, 

otherwise 

(35) 


The  Fourier  transform  of  f(x)  is: 


/(/?)  =  e 


2sm(/?^-)       f  sin(/?-^-)  -J2- 


/4 


0 


/4 


+ 


2A/ 


|sin(^-f-)-sin(^-^)j 


sin(jS -*-)}•     (36) 


and, 


e  =  ^-  +  A/(5-Q 


(37) 


The  characteristic  impedance  is: 


Zn  = 


vv  CC0 


(38) 


C 

^    "    Cn 


(39) 


where  v  =  3  *  10?  m  s,  and  C0  is  calculated  by  setting  zT  —  1. 

The  variational  expression  with  the  elementary  charge  distribution  function  f(x)  = 
|x|,  is  found  to  be  the  one  that  maximizes  the  value  of  capacitance  C,  thus  giving  the 
closest  value  to  the  exact  result  for  the  capacitance.  In  the  simple  distributions  chosen, 
M  is  the  variational  parameter,  and  maximum  capacitance  C  is  found  by  reducing  \I  to 
a  low  value  (see  page  18).  Figure  11  on  page  16  and  Figure  12  on  page  17  show  the 
calculated  results  of  impedance  and  effective  dielectric  constant  of  CBCPW  using  vari- 
ational method. 
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Figure  11.        Calculated    characteristic    impedance    of    CBCPW    by    variational 
method:    (e,  =  3.52,  H  =  1.5  mm) 
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VARIATIONAL  METHOD 
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Figure  12.       Calculated   effective   dielectric    constant   of   CBCPW   by  variational 
method:     (c,  =  3.52,  H  =  1.5  mm) 
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C.     COMPARISON  OF  RESULTS 

A  coplanar  waveguide  consists  of  a  strip  of  thin  metallic  film  on  the  surface  of  a 
dielectric  slab  with  two  ground  electrodes  running  adjacent  and  parallel  to  the  strip. 
Two  methods,  the  boundary  point  matching  method  and  the  variational  method,  were 
used  to  calculate  the  effective  dielectric  constant  (Ere//)  and  impedance  (Z0)  as  a  function 
of  the  line  dimensions. 

The  calculated  impedance  and  effective  dielectric  constant  of  conductor  backed 
coplanar  waveguide  are  shown  in  Figures  6,1,  11  and  12.  To  compare  the  results  of  the 
boundary  point  matching  method  and  variational  method,  several  graphs  are  plotted 
(Figure  13  on  page  20  through  Figure  20  on  page  27)  for  various  line  parameters. 

In  the  case  of  the  Rowe  and  Lao  boundary  point  matching  method  calculation,  the 
impedance  and  effective  dielectric  constant  are  computed  by  choosing  the  number  of 
matching-points,  X,  equal  to  100.  The  accuracy  of  this  computation  increases  with  N, 
but  the  computation  time  of  the  N  by  X  matrix  increases  approximately  as  A'4.  A  value 
of  X  =  100  is  chosen  in  the  present  work,  for  a  moderate  computation  time  and  accept- 
able accuracy  of  the  results.  Choice  of  larger  X  leads  to  roundoff  errors  and  excessive 
computation  time  in  the  matrix  calculation. 

The  characteristic  impedance  increases  as  the  center  conductor  width  decreases  and 
decreases  as  the  slot  widths  decrease.  As  shown  by  the  Rowe  and  Lao  boundary  point 
matching  method,  when  the  aspect  ratio  AH  >  5  the  slot  width  ratio  has  small  effect 
on  the  impedance  values.   This  behavior  is  shown  in  Figure  6  on  page  9. 

For  a  large  aspect  ratio  AH  versus  the  ratio  of  B/A,  the  impedance  (Z0)  does  not 
change  much,  since  both  the  slot  coupling  effect  and  the  effect  of  the  ground  backing 
are  small.  This  shows  that  characteristic  impedance  (Z0)  reflects  the  coupling  effect  be- 
tween the  two  slots.  Therefore,  by  adjusting  both  slot  widths  and  center  conductor 
width,  independent  control  of  the  characteristic  impedance  can  be  obtained.  For  the 
Rowe  and  Lao  boundary  point  matching  method,  the  effective  dielectric  constant  in- 
creases as  the  center  strip  width  increases  and  decreases  as  the  slot  widths  decrease. 

It  is  interesting  to  note  that  the  effective  dielectric  constant  treff  of  CBCPW  ap- 
proaches microstrip  behavior  in  its  dependence  on  the  linewidth  ratio  A/H,  especially  for 

wide  slots  (B/A  =  4.0  in  Figure  7  on  page  10).    For  microstrip,  as  microstrip  line  ratio 

£r  +  1 
A/H  goes  to  zero,  its  value  of  treff  approaches  — r—z —  =  2.26,  in  this  case,  and  as  A/H 

goes  to  oo,  then  ereff  goes  to  £,.    Similar  behavior  is  also  noted  for  CBCPW  in  Figure  7 

on  page  10. 
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In  case  of  Yamashita's  variational  method,  the  impedance  and  the  effective  dielectric 
constant  are  calculated  using  a  value  of  M  =  0.0-4.  The  ratio  of  the  charge  density 
magnitude  value  between  the  center  strip  and  outer  strip  is  thus  in  the  range  20:1  -  50:1. 
The  impedance  and  the  effective  dielectric  constant  are  compared  with  the  Rowe  and 
Lao  results  by  calculation  using  the  same  parameters. 

The  results  from  the  variational  method  are  in  good  agreement  with  those  of  the 
Rowe  &  Lao  method.  For  impedance,  Figure  13  on  page  20  through  Figure  16  on  page 
23  shows  that  less  than  6%  difference  in  Z„  was  observed  for  most  cases,  and  also 
Figure  17  on  page  24  through  Figure  20  on  page  27  shows  less  than  5%  deviation  in 
the  effective  dielectric  constant  values.  The  outputs  of  the  both  methods  are  closer  as 
the  slot  widths  increase.  The  characteristics  approach  those  of  microstrip  as  the  slot 
widths  increase  and  the  ground  planes  have  less  effect.  For  example,  microstrip  lines 
on  this  substrate  having  linewidth  ratios  AH  equal  to  1  and  10  would  have  character- 
istic impedances  Z0  equal  to  78.7  ohms  and  16.6  ohms,  respectively,  and  effective 
dielectric  constant  treff  of  2.61  and  3.13,  respectively.  These  are  comparable  to  the 
wide-slot  values  (B  A  =  4.0)  for  CBCPW  seen  in  Figures  16  and  20.  The  computation 
time  on  the  IBM  3278  for  one  data  point  is: 

•  Boundary  point  matching  method  =  50  sec. 

•  Variational  method  =   10  sec. 
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Figure  13.  Computed  results  for  characteristic  impedance  versus  strip  width  with 
line  dimensions  as  parameters:  (  B/A=  1.5,  C/B  =  2.0,  er  =  3.52,  H  = 
1.5  mm ) 


20 


COMPARISON  OF  METHODS 

IMPEDANCE 

o 

CO 

—•—  ROWE  AND  LAO 
— * —  VARIAT10NAL 

B/A=2.0,  C/B=2.0 

o 

B 

D 

N 

? 

o 

I                 i 

I                  I                  I                 I                  I                  I                  I 

i 

2 

4                                6                                B 

A/H 

10 

Figure  14.  Computed  results  for  characteristic  impedance  versus  strip  width  with 
line  dimensions  as  parameters:  (  B,'A=2.0,  C,B  =  2.0,  zr  =  3.52,  H  = 
1.5  mm  ) 
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Figure  15.  Computed  results  for  characteristic  impedance  versus  strip  width  with 
line  dimensions  as  parameters:  (  B,A  =  3.0,  C/B  =  2.0,  cr  =  3.52,  H  = 
1.5  mm  ) 
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Figure  16.  Computed  results  for  characteristic  impedance  versus  strip  width  with 
line  dimensions  as  parameters:  (  B/A  =  4.0,  C,B  =  2.0,  cr  =  3.52,  H  = 
1.5  mm  ) 
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Figure  17.       Computed  results  for  dielectric  constant  versus  strip  width  for  various 
line  parameters:     (  B.'A=  1.5,  C/B  =  2.0,  er  =  3.52,  H  =  1.5  mm  ) 
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Figure  18.      Computed  results  for  dielectric  constant  versus  strip  width  for  various 
line  parameters:     (  B,'A  =  2.0,  C/B  =  2.0,  tr  =  3.52,  H  =  1.5  mm  ) 
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Figure   19.       Computed  results  for  dielectric  constant  versus  strip  width  for  various 
line  parameters:     (  B/A=3.0,  C/B  =  2.0,  er  =  3.52,  H  =  1.5  mm  ) 
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Figure  20.      Computed  results  for  dielectric  constant  versus  strip  width  for  various 
line  parameters:     (  B,A  =  4.0,  C/B=2.0,  z,  =  3.52,  H  =  1.5  mm  ) 
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III.     CONSTRUCTION  OF  COPLANAR  LINES  AND  TRANSITIONS 

A.     ADHESIVE  FOIL  METHOD 

In  this  section,  the  adhesive  foil  method  for  circuit  fabrication  is  described.  For  the 
design  of  a  coplanar  line  it  is  essential  to  have  accurate  knowledge  of  the  dielectric  con- 
stant of  the  substrate.  An  illustration  of  the  use  of  adhesive  foil  in  planar  circuitry  is  in 
the  measurement  made  to  determine  the  substrate  er  by  means  of  a  microstrip  resonator 
circuit. 

To  measure  the  dielectric  constant  of  substrate  used  in  the  present  work,  a  quarter 
wave  length  stub  of  conventional  microstrip  is  used  as  shown  in  Figure  21. 


One  quarter  wave  length  microstrip 

Substrate  height  =  1.5  mm 

A  /  H  (  Line  width  ratio  )  =  2.28 


L  = 


ms 


/  4 


Figure  21.      Quarter  wave  length  stub  for  substrate  dielectric  constant  measurement 

The  line  pattern  is  fastened  to  the  substrate  using  adhesive  copper  foil.    The  value  of 
S2l  or  S[2,  transmission  coefficient,  is  measured  with  respect  to  frequency  using  the  net- 
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work  analyzer  to  find  the  resonant  frequency,  which  occurs  at  maximum  dB  transmission 
loss  (see  Figure  22  on  page  29). 

S2i 

A 


\  / 


->- 


f0  Frequency 


Figure  22.      DB  Magnitude  plot  for  S2l 


The  following  equations  are  used  to  find  the  dielectric  constant  of  the  substrate: 


/  = 


'■ms 


{at  f0) 


(40) 


C 


(41) 


therefore, 


-  (      c      \2 

JO'-ms 
=  (  — )2 

{  Wo  } 


3x  10s 


4x0.03x  1.5  x  101 

=  2.8 


(42) 


where  zrtm  is  the  microstrip  value  of  effective  dielectric  constant  on  the  substrate  material 
employed.  Using  standard  microstrip  relations,  the  substrate  relative  dielectric  constant 
was  found  to  be:   e,  =  3.52. 
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And  next,  with  this  dielectric  constant  (er),  the  CBCPW  parameters  may  be  com- 
puted by  using  the  Rowe  and  Lao  method  which  is  presented  in  the  previous  chapter, 
and  the  coplanar  line  may  be  constructed  on  the  same  substrate  material. 

The  final  steps  for  fabrication  of  CBCPW  are  the  attachment  of  the  line  strips  to  a 
transition  to  coaxial  line  and  solderins  the  outer  conductors  to  around  with  good  con- 
nections,  to  reduce  the  insertion  loss.  In  this  adhesive  foil  method,  the  disadvantages 
are  lack  of  precision  and  difficult  assembly,  since  cutting  and  attaching  is  done  manually. 

B.     PHOTOLITHOGRAPHICAL  METHOD 

To  overcome  the  inaccuracy  of  the  adhesive  foil  method,  a  photolithographical 
method  for  CBCPW  design  is  used.  This  method  employs  a  CAD  program, 
photolithography  and  an  etching  process  to  form  the  line  pattern. 

Real-size  design  is  possible  on  the  computer  by  using  the  available  CAD  software 
tool  with  CBCPW  parameters  from  the  Rowe  and  Lao  calculation.  A  real-size  negative 
film  is  then  made  at  the  photo  shop  (This  work,  was  done  at  the  XPS  Educational  VIedia 
Department).  The  exposure  of  the  CBCPW  pattern  is  done  with  this  negative  film. 
Figure  23  on  page  31  shows  drawn  real-size  line  pattern  and  negative. 


Figure  23.      Drawn  line  pattern  and  negative 


The  circuit  board  etching  process  is: 
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1.  Expose  photosensitive  emulsion  through  the  prepared  mask. 

2.  Allow  board  to  run  through  the  etcher  to  remove  all  unwanted  copper. 

3.  After  ensuring  by  visual  inspection  that  all  unwanted  copper  has  been  removed, 
remove  negative  film  from  board  using  acetone. 

The  advantage  of  this  method  is  that  greater  precision  is  achieved  than  with  the 
adhesive  foil  method. 

C.     TRANSITION  TO  COAXIAL  LINE 

It  is  necessary-  to  make  a  transition  from  the  CBCPW  to  a  conventional  transmission 
medium.  An  important  factor  in  making  such  transitions  is  an  impedance  match  and  a 
field  match.  The  field  match,  which  is  the  prerequisite  to  a  proper  impedance  match,  is 
provided  by  the  mechanical  transition  and  the  impedance  match  is  provided  by  the 
electrical  design.  In  this  section,  the  mechanical  transition  of  the  CBCPW  to  coaxial  line 
will  be  described. 

The  transition  from  a  microstrip  line  to  coaxial  line  is  obtained  by  connecting  the 
center  pin  of  the  coaxial  connector  to  the  microstrip.  In  contrast  to  the  microstrip  line, 
a  conductor  backed  coplanar  waveguide  has  ground  conductors  and  a  center  conductor 
(signal  line)  on  the  one  side  of  the  substrate,  and  it  also  has  a  ground  plane  on  the  other 
side  of  the  substrate.  Therefore,  the  requirement  of  CBCPW  transition  to  coaxial  line 
is  different  than  the  microstrip  case. 

The  CBCPW  transition  to  coaxial  line  is  obtained  by  extending  the  center  pin  of  the 
coaxial  connector  to  the  center  conductor  of  CBCPW  to  make  signal  contact,  and  by 
connecting  the  ground  flange  of  the  coaxial  connector  to  outer  conductor  and  back- 
plane of  CBCPW.  Therefore,  top  and  bottom  ground  plane  connectors  are  required  to 
make  a  transition  from  a  CBCPW  to  coaxial  line. 

The  objective  of  the  coaxial  line  to  CBCPW  transition  is  to  obtain  an  impedance 
match  and  very  low  insertion  loss.  When  attaching  the  coaxial  connector  to  the 
CBCPW,  it  is  essential  to  have  a  good  contact  which  is  as  direct  (short)  as  possible  to 
reduce  insertion  loss.  This  requires  careful  mechanical  design  to  conform  to  the  geom- 
etry of  the  CBCPW. 

A  coupling  unit  was  constructed  to  make  the  transition  from  CBCPW  to  N-type 
coaxial  line.  In  addition  to  providing  the  connection  from  the  coaxial  center  conductor 
to  the  CBCPW  center  conductor,  a  pair  of  tab  connectors  were  provided  at  the  top  of 
the  substrate  so  that  soldered  ground  connections  from  the  fixture  to  the  two  coplanar 
ground  lines  could  be  made. 
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IV.     MEASUREMENT  OF  COPLANAR  LINE  PARAMETERS 

A.     RESONATOR  METHOD  FOR  PHASE  VELOCITY  AND  EFFECTIVE 
DIELECTRIC  CONSTANT 

In  this  section,  two  methods  for  calculating  phase  velocity  and  effective  dielectric 
constant  were  explored.  The  first  method  used  a  series  open-ended  resonator;  the  second 
method  used  a  shunt  resonator. 

1.     Series  Open-ended  Resonator 

A  series-connected  resonator  is  constructed  using  a  two  port  section  of  CBCPW 
connected  via  coaxial  connectors.  As  shown  in  Figure  24  on  page  34,  gaps  in  the 
CBCPW  center  strip  are  separated  by  a  half  wavelength  interval  of  coplanar  line.  At  this 
length  the  open-ended  section  will  act  as  a  series-connected  resonator.  A  gap  length  of 
a  few  tenths  of  a  line  width  provides  sufficient  coupling,  with  minimum  resonator  load- 
ing. 

The  equivalent  circuit  model  for  a  series  open-ended  resonator  is  shown  in  Fig- 
ure 24  on  page  34  as  designed  for  use  with  the  Touchstone  CAD  program.  In  this  way, 
we  can  compare  the  Touchstone  model  for  CBCPW  with  the  results  of  laboratory 
measurement.  Measurements  were  made  to  find  the  resonant  frequency  by  measuring 
the  transmission  matrix  element  S2i  in  dB  with  respect  to  frequency.  The  resonant  fre- 
quency occurs  at  minimum  dB  transmission  loss  as  shown  in  Figure  25  on  page  35. 

The  following  equations  are  used  to  find  the  effective  dielectric  constant  of 
CBCPW  by  using  the  resonant  frequency  from  a  series  open-ended  resonator. 

1  =  -^-        (atf0)  (43) 

vcpw      JO^cpw 
C 

\'Cre/f  (44) 

jgjg 

\ereff 
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a  :    Series  open-ended  resonator 
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b  :    Equivalent  circuit       (  C  =  0.005  pF 


Figure  24.      Series  open-ended  resonator  and  equivalent  circuit 
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where  Xcpw  and  vcpw  are  the  wavelength  and  phase  velocity  of  the  waves  on  the  CBCPW 
and  /„  is  the  free  space  wavelength. 


EEsof  -  Touchstone  -  Tue  Nov  01  23:  00:  05  19B8  -  JE0NG11 

□   DB[S21] 
CBCPW 


0.000 


-60.  00 


-120.  0 


1.000 


-e b b- 


FHEQ-GHZ 


1.500 


Figure  25.      Plot  for  dB  magnitude  Slx  versus  frequency 

The  value  of/  used  in  Touchstone  CAD  program  is  92.8  mm  at  a  frequency  of 
1.0  GHz.  Thus  the  calculated  value  of  the  effective  dielectric  constant  at  A/H  =  2.0, 
B  A  =  2.0,  C/B  =  1.5  is: 


Zreff-(i^)2 


3x  10a 


2  x  0.0928  x  1.0  x  109 
=  2.61 


(46) 
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2.     Shunt  Resonator 

Two  ways  to  calculate  the  resonant  frequency  by  using  a  Touchstone  CAD 
program  for  a  CBCPW  resonator  are  explored.  One  is  a  one-port  open-ended  resonator; 
the  other  is  a  shunt  short-circuited  resonator. 
a.     One-port  Open-ended  Resonator 

In  Figure  26  on  page  36,  the  equivalent  circuit  is  shown  for  an  open-ended 
stub  resonator  for  the  input  impedance  calculation  by  the  Touchstone  CAD  program. 
The  length  of  the  strip  is  a  quarter  wavelength.  The  open-ended  stub  input  impedance 
is: 


=         ^CPw/ 


Figure  26.      Equivalent  circuit  for  calculation  of  angle  of  impedance 


Zin  =  -jZ0  cot(/?/) 
=  -jZ0  cot(  — - ) 


(47) 


here  Zm  changes  from  a  negative  to  a  positive  value  at  resonance,  that  is,  the  angle  of 
Zm  changes  from  — '—  radians  to  +  -z-  radians.  Therefore,  the  difference  of  angle  of  in- 
put impedance  is  n  radians.  This  behavior  is  plotted  in  Figure  28  on  page  39  using  the 
Touchstone  CAD  program. 

The  equations  to  calculate  the  effective  dielectric  constant  for  the  open- 
ended  resonator  are  the  same  as  for  the  series  open-ended  resonator  case  except  that 
here,  /=/e,w/4. 

Therefore,  the  Final  equation  for  calculation  of  the  effective  dielectric  con- 
stant is: 
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=  (A)2  (48) 

1  4/  ; 

=  (  —  ? 

where  c  is  the  velocity  of  light  and^  is  the  resonant  frequency. 

The  value  of/  used  in  the  Touchstone  CAD  program  is  15.45  mm  to  obtain 
a  -j-  resonator  at  a  frequency  of  3  GHz,  based  on  the  result  of  the  series  open-ended 
resonator  case.  Thus  the  calculated  value  of  the  effective  dielectric  constant  at  A/H  = 
2.0,  B/A  =  2.0,  C/B  =  1.5  is: 

:reff: 

-"■JO 

c  0.01545  x  3.0  x  109  ' 
=  2.61 

b.     Shunt  Short-circuited  Resonator 

As  shown  in  Figure  27  on  page  38,  a  short-terminated  section  of  CBCPW 
having  a  length  of  one-half  coplanar-guide  wavelength  is  connected  in  shunt  with  the 
through  line.  This  is  the  shunt  shorted  resonator  which  models  the  experimental  circuit 
shown  in  Figure  23  on  page  31.  A  transmission  minimum  occurs  at  resonance.  The 
resonant  frequency  of  this  shunt  resonator  can  also  be  calculated  using  the  Touchstone 
CAD  program,  with  use  of  an  equivalent  circuit  model.  (See  Figure  27  on  page  38.)  The 
input  impedance  of  shunt  shorted  stub  is: 

Zin  =  -jZ0  tan(/?/) 

■7  t     (  2nt  s  (50) 

-  -jZ0  tan(  — ) 

When  /  is  equal  to  -f-,  the  input  impedance  goes  to  zero.  Therefore,  a  minimum  value 
of  S21  can  be  obtained  from  Touchstone  CAD  program  at  the  resonant  frequency.  The 
resonant  frequency  at  the  minimum  value  of  S2l  is  then  used  for  calculation  of  the  ef- 
fective dielectric  constant  of  the  CBCPW. 
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a  :    Shunt  shorted  resonator 
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b  :  Equivalent  circuit 


nj 


Figure  27.      Shunt  shorted  resonator  and  equivalent  circuit 

The  dB  magnitude  output  of  the  line  section  with  the  shunt  shorted 
resonator  is  plotted  in  Figure  28  on  page  39  and  compared  with  the  plot  of  the  angle 
of  input  impedance  of  the  open-ended  stub.  The  equations  to  calculate  the  effective 
dielectric  constant  are  the  same  in  form  as  the  series  open-ended  resonator,  since  these 
are  both  half-wavelensth  resonators. 
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Figure  28.      Plot  for  dB  magnitude  S2l  with  shunt  shorted  stub,  and  computed  angle 


of  impedance  of  -j-  open  stub  at  resonance 
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Table  1  on  page  40  contains  CBCPW  effective  dielectric  constants  as  cal- 
culated from  the  results  of  the  two  types  of  measurement  carried  out:  those  with  the 
series  -f-  resonator  and  the  shunt  -f-  resonator.  The  CBCPW  used  in  these  measure- 
ments had  dimensions:  A  =  3  mm,  B  =  6  mm,  C  =  9  mm. 

Table   1.     COMPARISON  OF  THE  EFFECTIVE  DIELECTRIC  CONSTANT 


Item 

Effective  dielectric  constant  (f0) 

Touchstone 

Measured 

Series  open-ended  stub 
(1  =  92.8  mm) 

2.61  (1  GHz) 

2.61  (1  GHz) 

Shunt  shorted  stub 
(1  =  30.9  mm) 

2.61  (3  GHz) 

2.68  (2.96  GHz) 

B.     LEAKAGE  WAVE  MEASUREMENTS  ON  COPLANAR  LINE  AND 
MICROSTRIP 

Conductor  backed  coplanar  waveguide  (CBCPW)  transmission  line  is  compared  to 
microstrip  in  terms  of  relative  extent  of  the  ambient  fields  around  the  lines.  An  angle 
section  is  used  so  that  the  fields  both  around  a  corner  and  straight  sections  may  be  ob- 
served. The  angle  section  is  expected  to  be  a  worst  case  due  to  the  greater  magnetic- 
energy  storage  in  the  neighborhood  of  the  corner. 

1.     Measurement 
The  steps  are: 

1.  Choose  the  dimensions  of  CBCPW  for  50  Q  impedance  by  using  the  Rowe  and  Lao 
method  (substrate  height  is  H  =1.5  mm,  dielectric  constant  g,  =  3.52). 

2.  Design  an  angle  stub  (See  Figure  29  on  page  41)  to  measure  the  fields  near  the 
conductor  at  the  angle  point. 

3.  Make  a  small  loop  antenna  to  detect  the  leakage  wave  (diameter  =   10  mm)  (see 
Figure  30  on  page  41). 
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Figure  29.      Line  angle  for  conductor  ambient  field  measurement 
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Figure  30.      Block  diagram  for  conductor  ambient  field  measurement 
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4.  Move  this  probe  (loop  antenna)  in  a  straight  line  from  the  middle  point  of  the 
center  conductor  of  the  angle  portion  to  outside  of  the  line. 

5.  Read  the  oscilloscope  output  and  plot  the  results  to  compare  the  outputs  of 
microstrip  and  CBCPW. 

Although  as  is  true  in  all  such  field-probing  measurements,  the  electromagnetic 
fields  are  perturbed  by  the  presence  of  the  probe,  and  no  field  details  smaller  than  the 
size  of  the  probe  can  be  observed,  a  useful  comparison  of  the  two  cases  (coplanar  line 
and  microstrip)  observed  can  be  obtained. 
2.     Results 

Figure  31  on  page  43  shows  the  probe  measurement  plotted  against  the  distance 
for  microstrip  and  CBCPW  on  1.5  mm  thick  substrate  with  a  relative  permittivity  of  3.52 
at  3  GHz.  A  thin  copper  conductor,  which  is  designed  by  the  photolithographical 
method,  is  chosen. 

For  a  given  substrate  thickness  and  dielectric  constant,  CBCPW  has  a  lower 
ambient  field  strength  than  microstrip.  The  spurious  coupling  between  adjacent 
coplanar  lines  is  therefore  expected  to  be  significantly  lower  than  with  microstrip. 
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Figure  31.      Conductor  ambient  field  distribution 
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V.     SUMMARY  AND  CONCLUSION 

In  the  foregoing  investigation,  conductor  backed  coplanar  lines  on  finite-height 
substrate  have  been  investigated.  (In  this  present  work,  substrate  height  H  =  1.5  mm, 
dielectric  constant  zr  =  3.52.)  This  is  a  typical  form  for  use  in  microwave  integrated 
circuitry. 

The  research  has  included: 

1.  Calculation  of  the  effective  dielectric  constant  and  characteristic  impedance  of 
conductor  backed  coplanar  waveguide. 

2.  Construction  of  coplanar  lines  and  a  transition  from  CBCPW  to  a  coaxial  connec- 
tor. 

3.  Measurements  of  CBCPW  parameters  (phase  velocity  and  effective  dielectric  con- 
stant) on  practical  line  sections. 

Two  computational  methods,  the  boundary  point  matching  method  and  the  vari- 
ational method,  were  used  to  calculate  the  effective  dielectric  constant  and  impedance 
for  various  line  dimensions.  The  impedances  and  effective  dielectric  constant  obtained 
from  the  variational  method  are  compared  with  values  obtained  from  the  boundary 
point  matching  method  as  computed  for  similar  line  parameters.  The  results  of  both 
methods  are  in  good  agreement.  For  impedance  (Z0),  Figure  13  on  page  20  through 
Figure  16  on  page  23  show  that  less  than  6  %  difference  was  observed  for  most  cases, 
and  less  than  5  %  difference  in  the  dielectric  constant  is  shown  in  Figure  17  on  page  24 
through  Figure  20  on  page  27. 

For  a  large  aspect  ratio  AH,  impedance  shows  little  change  versus  the  ratio  B  A. 
since  both  slot  coupling  effects  and  the  effects  of  fringing  fields  are  relatively  small  with 
large  A.  H  ratio.  A  wider  gap  causes  a  larger  impedance,  and  the  microstrip  character- 
istic behavior  will  be  approached  when  the  gap  goes  to  infinity  or  zero.  Therefore,  by 
adjusting  both  slot  widths  and  center  strip  width,  independent  control  of  the  character- 
istic impedance  can  be  obtained. 

Figure  32  on  page  46  and  Figure  33  on  page  47  show  that  the  impedance  and  the 
effective  dielectric  constant  of  CBCPW  are  changed  appreciably  with  increasing 
dielectric  constant  of  the  substrate.  The  choice  of  proper  dielectric  constant  of  substrate 
is  therefore  very  important.  A  method  for  calculating  this  dielectric  constant  of  the 
substrate  material  is  detailed  in  Chapter  III  (A). 
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As  was  mentioned  in  Chapter  I  (Introduction),  CBCPW  has  a  smaller  extent  of 
ambient  fields  than  does  microstrip.  The  spurious  coupling  between  coplanar  lines  and 
adjacent  circuitry  is  thus  lower  than  with  microstrip.  This  behavior  is  demonstrated  in 
Chapter  IV,  through  leakage  field  measurements  on  coplanar  line  and  microstrip. 

In  Chapter  III,  two  methods  for  fabrication  of  CBCPW  are  considered,  and  we  infer 
that  the  photolithographical  method  is  better  than  the  adhesive  foil  method  due  to  its 
inherently  greater  accuracy.  Also  the  three  special  requirements  of  CBCPW  transition 
to  coaxial  line  are  described. 

These  three  requirements  are: 

1.  Top  and  bottom  ground  plane  connector. 

2.  Connection  lead  length  as  short  as  possible. 

3.  Smooth  mechanical  transition  between  physical  components  to  avoid  parasitics. 

Phase  velocity  and  effective  dielectric  constant  of  CBCPW  can  be  calculated  by  us- 
ing two  kinds  of  resonator:  series  open-ended  and  shunt  shorted,  as  described  in  Chap- 
ter IV.  This  consideration  gives  the  following  points.  The  series  resonator  must  be 
lightly  coupled  in  order  to  avoid  pulling  of  the  resonant  frequency.  In  the  present  work 
the  parasitic  capacitance  to  ground  of  the  open  ends  was  unaccounted  for,  but  these 
values,  as  inferred  from  the  corresponding  microstrip  case,  are  typically  small. 

The  availability  of  design  data  for  conductor  backed  coplanar  line  will  facilitate  the 
use  of  this  transmission  format  in  microwave  integrated  circuit  design.  It  is  anticipated 
that  the  results  of  this  study  will  promote  a  better  understanding  of  coplanar  waveguide, 
and  promote  its  use  in  practical  microwave  circuits. 
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ROWE  AND  LAO  METHOD 
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Figure  32.      Characteristic  impedance  versus  strip  width  with  dielectric  constant  as  a 
parameter:     (  B'A=2.0,  CB  =  2.0,  H  =  1.5  mm  ) 
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ROWE  AND  LAO  METHOD 
EFFECTIVE  DIELECTRIC  CONSTANT 
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Figure  33.      Effective  dielectric  constant  versus  strip  width  with  substrate  dielectric 
constant  as  a  parameter:    (  B/A  =  2.0,  C/B  =  2.0,  H  =  1.5  mm  ) 
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APPENDIX  A.     FORTRAN  PROGRAM  TO  CALCULATE  THE 

IMPEDANCE  AND  EFFECTIVE  DIELECTRIC  CONSTANT  OF  CBCPW 

USING  BOUNDARY  POINT  MATCHING  METHOD 

C$NOEXT 

r  *******     PURPOSE  **i 


C     TO  CALCULATE  THE  IMPEDANCE  AND  EFFECTIVE  DIELECTRIC  CONSTANT 
C     OF  CBCPW  BY  USING  BOUNDARY  POINT  MATCHING  METHOD 


C  *******     VARIABLE  DEFINITIONS  ******* 

C  H     SUBSTRATE  HEIGHT 

C  ER  =  DIELECTRIC  CONSTANT  OF  SUBSTATE 

C  ERE  =  EFFECTIVE  DIELECTRIC  CONSTANT  OF  CBCPW 

C  CI  =  CAPACITANCE  FOR  ER  =  3.52 

C  CO  =  CAPACITANCE  FOR  ER  =  1.0 

C  ZO  =  IMPEDANCE  OF  CBCPW 

C  V     VELOCITYOF  LIGHT 

C  AA  =  RATIO  OF  A/H 

C  BB  =  RATIO  OF  B/A 

C  CC  =  RATIO  OF  C/B 


REAL*8  H,ER,K,M,A,B,C,G,L,C0,C1,Z0,ERE,V,AA,BB,CC 

V=3.D+08 

0PEN(UNIT=7,FILE='RLM  DATA') 


10  PRINT* 
READ(* 
PRINT* 
READ(* 
PRINT* 
READ(* 
PRINT* 
READ(* 
PRINT* 
READ(* 
PRINT* 
READ(* 


'HOW  MANY  POINTS  FOR  CALCULATING  THE  IMPEDENCE  (ZO)?' 
■)N 

'THICKNESS  OF  SUBSTRATE  (H(M))?' 

)H 

' ER=? ' 
■)ER 

'A/H=?' 
•)AA 

'B/A=?' 
•)BB 

'C/B=?' 
■)CC 


A=AA*H/2.  D00 
B=BB*A 
C=CC*B 
L=2. D00*C 

CALL  LAO(N,H,l.D00,A,B,C,L,C0) 

CALL  LA0(N,H,ER,A,B,C,L,C1) 
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ERE=C1/C0 

Z0=( 1. D00/V)/SQRT(C1*C0) 

WRITE(7,33)AA,BB,CC,Z0,ERE 
33  F0RMAT(3(1X,D8.4),2(1X,D9.5)) 

PRINT'S    'IF  YOU  WANT  TO  RUN  AGAIN, ENTER  (22)' 

READ(*,*)MM 

IF  (MM.  EQ.  22)  GOTO   10 

STOP 

END 


/""I  «J-.  mJm  «.'»  uim  m3m  >J*  «.'»  «J*  «J j  i.'*  «.'»  aJa  -J*  m3m  mjf  aJLa  <JLa  «J.  - '»  -.'-  «J-  -J-  — '-•  «*.»  »' -  J-  «.'-  «J-  iJLa  aaTa  -*.•  J—  -.'.. 

SUBROUTINE   LA0(N,H,ER,A,B,C1,L,C) 


REAL*8  K(N) ,X(N) ,M(N,N) ,T(N,N+1) ,D(N) ,AA(N) 
#  ,H,ER,A,B,G,L,Z0,ERE,PHI,E0,C0,C,V,C1 

PHI=4.  D00*DATAN(  1.  D00) 

E0=lD-09/(36*PHI) 

C=0 

DO   100   1=1, N 

K(I)=(2*I-1)*PHT/(2*L) 
DO   200   J=1,N 

X(J)=J*L/(N+1) 

IF((X(  J).  GT.  A.  AND.  X( J).  LT.  B).  OR.  X(  J).  GT.  C1)THEN 

M(J,I)=K(I)*C0S(K(I)*X(J))*(1+ER*1/TANH(K(I)*H)) 
ELSEIF(  (X(  J).  GE.  B.  AND.  X(  J).  LE.  CI).  OR.  X(  J).  LE.  A)THEN 

M(J,I)=COS(K(I)*X(J)) 
ELSE 

PRINT- , ' SORRY,   WRONG  APPLY.    CHECK   INPUT  DATA  AGAIN!  * 
ENDIF 

T(J,I)=M(J,I) 
IF(X(J).GT.  A)D(J)=0 
IF(X(J).LE.  A)D(J)  =  1 
200  CONTINUE 

T(I,N+1)=D(I) 
100        CONTINUE 


CALL  GAUSS(T,N,N+1,AA) 

DO  300   1=1, N 

C=AA(I)*SIN(K(I)*(A+B)/2)*(1+ER*1/TANH(K(I)*H))+C 

300        CONTINUE 

C=-2. D00*E0*C 

RETURN 
END 
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»'-  «  -'-  ^'-  ~*-  »-  »'-  »'-  J-  J-  J-  -•-  J-*'-  .a.  »*-  J-  J-  «;*  -;-  j,  j.  »•..  „*.  jl 

SUBROUTINE   GAUSS ( A,N,L,X) 

JhJUJUULJU^L ■  * ■  ■  T  ■   f  i -.'..i. -J- ^. .'.. .J^  a.  -A. »f, .J-  J^J. .J. .j. .j, .j,. 


C  SOLUTION  BY  GAUSSIAN  ELIMINATION   

REAL*8  A(N,N+1),X(N) 

NN=N-1 

DO  99   K=1,NN 

CALL  FIND(A,N,N+1,K,J) 
CALL  CHANGE(A,N,N+1,K,J) 
KK=K+1 

DO  88   I=KK,N 

CALL  ROWMUL(A,N,N+l,K,I,-A(I,K)/A(K,K)) 
88  CONTINUE 

99  CONTINUE 

CALL  BKSUB(A,N,N+1,X) 

RETURN 
END 


SUBROUTINE  FIND(A,N,L,K, J) 

JL.  A  J-  J-  JUJUJUJa  ^'-  .J-  .A.  .A..A.. I-  •.*..  **--.■->•-  JL  >*--.•-  JL  JUJL  J-  .J. 


C--   FIND  THE  LARGE  VALUE (TO  AVOID  THE  FLOATING  POINT  DIVIDED  ZERO) 

REAL*8  A(N,N+1) 

J=K 

KK=K+1 

DO   100   I=KK,N 

IF   (ABS(A(J,K)).LT.ABS(A(I,K)))    J=I 
100     CONTINUE 

RETURN 
END 


SUBROUTINE  CHANGE ( A,N,L,K, J) 


C--   CHANGE  THE  ROW   (TO  AVOID  THE  FLOATING  POINT  DIVIDED   ZERO   ) 
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100 


c 


REAL*8  A(N,N+1) 
IF   (J.  EQ.K) RETURN 


DO   100   1=1 
T=A(K, 
A(K,I> 
A(J,I)' 

CONTINUE 

,L 

I) 

=A(J: 

=T 

,D 

RETURN 
END 

SUBROUTINE  ROWMUL(A,N,L,K, J,D) 


C  ROW  MULTIPLICATION   -- 

REAL* 8  A(N,N+1),D 

DO   100   1=1, L 

A(J,I)=A(J,I)+D*A(K,I) 
100       CONTINUE 

RETURN 
END 


p        ju  j-  ju  -»-  j-  j-  j-  j-  ju  ju  ju  ju  j-  ju  J-  J-  ju  j-  ju  ju  j-  j -  ju  J-  ju 

SUBROUTINE   BKSUB( A,N,L,X) 


U  J-  .'-  JU  JU  JU  J-  JU  JU  J-  JU  JU  J-  J.  J.  J*  J.  JU  JU 


C  FIND  THE  RESULT  OF  X(K)    - 

C 

REAL*8  A(N,N+1),X(N) 

X(N)=A(N,N+1)/A(N,N) 

NN=N-1 

DO   100  J=1,NN 
K=N-J 

C  CALCULATING  SUM  FOR  X(K) 

SUM=0. 0 
NNN=N-K 

DO  200  JJ=1,NNN 
M=K+JJ 

SUM=SUM+A(K,M)*X(M) 
200  CONTINUE 

C  SOLVING  FOR  X(K) 
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X(K)=(A(K,N+1)-SUM)/A(K,K) 
100   CONTINUE 

RETURN 
END 
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APPENDIX  B.     FORTRAN  PROGRAM  TO  CALCULATE  THE 
IMPEDANCE  AND  EFFECTIVE  DIELECTRIC  CONSTANT  OF  CBCPW 

USING  VARIATIONAL  METHOD 

C$NOEXT 

q  *******  PURPOSE  ******* 

C     TO  CALCULATE  THE  IMPEDANCE  AND  EFFECTIVE  DIELECTRIC  CONSTANT 
C     OF  CBCPW  BY  USING  VARIATONAL  METHOD. 

C     *******  VARIABLE  DEFINITIONS  ******* 

C  H     SUBSTRATE  HEIGHT 

C  ER  =  DIELECTRIC  CONSTANT  OF  SUBSTATE 

C  ERE  =  EFFECTIVE  DIELECTRIC  CONSTANT  OF  CBCPW 

C  CD  =  CAPACITANCE  FOR  ER  =  3.  52 

C  CO  =  CAPACITANCE  FOR  ER  =  1.0 

C  ZO  =  IMPEDANCE  OF  CBCPW 

REAL  Q0,Q,A,B,C,ER,ERE,PI,H,M1,CD,C0,Z0 
EXTERNAL  F 

COMMON  H,A,B,C,Q0.Q,PI,M1.ER 
0PEN(UNIT=2,FILE=fVM  DATAf) 

DATA  AERR,RERR,Al,Bl/l.E-06,0,0.  lE-30,3. 4E03/ 

PI=4.*ATAN(1.  ) 
Q0=1.0 
Ml=4. OE-02 
H=0.  0015 
ER=3.  52 

A=2. 1*H 
B=l.  0*A 
C=2. 0*B 

Q=Q0*A/2.+Ml*(B-C) 

S1=PI*8.  84*1.  E-12*Q**2 

S2=DCADRE(F,A1,B1,AERR,RERR, ERROR, IER) 

ER=1.  0 

S3=DCADRE(F,A1,B1,AERR,RERR, ERROR, IER) 

CD=S1/S2 
C0=S1/S3 
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ERE=CD/CO 

Z0=1.  /(3.*1.E08*SQRT(CD*C0)) 

WRITE(2,11)A/H,B/A,C/B3Z0,ERE 
11    F0RMAT(1X,5(F14.  8)) 
STOP 
END 


REAL  FUNCTION  F(BB) 

REAL  Y1,Y2,A,B,C,H,Q0,ER,BB,M1 
COMMON  H,A,B,C,Q0,Q,PI,M1,ER 

Yl=(Q*(4.*SIN(BB*A/(2. ))/(BB*A)-(4.*SIN(BB*A/(4.  ))/ 
#(BB*A))**2)+2*Ml/BB*(SIN(BB*B/2.  )-SIN(BB*C/2.  )))**2 

Y2=BB*( 1+ER/TANH( BB*H) ) 

F=Y1/Y2 

RETURN 
END 
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